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We discuss magnetoexcitons dispersion in single and coupled GaAs — {Ga, Al)As 
quantum wells using the Bethe-Salpeter (B-S) formalism. The B-S formalism in the case of 
quantum wells provides an equation for the exciton wave function which depends on two 
space variables plus the time variable, i.e. the B-S equation is 2 + 1-dimensional equation. 
We compare the results for magnetoexcitons dispersion, obtained in the LLL approxima- 
tions with the results calculated by solving the exact B-S equation. It is shown that the exact 
B-S equation has an extra term (B-S term) that does not exist in the LLL approximation. 
Within the framework of the variational method, we obtain that, (i) the ground-state energy 
of a heavy-hole magnetoexciton with a zero wave vector in GaAs — {Ga, Al)As quantum 
wells, calculated by means of the exact B-S equation, is very close to the ground-state en- 
ergy, obtained in the LLL approximation, (ii) in a strong perpendicular magnetic field the 
magnetoexciton dispersion (in-plane magnetoexciton mass) is determined mainly by the B- 
S term rather than the term that describes the electron-hole Coulomb interaction in the LLL 
approximation. 



1 Schrodinger equation for magnetoexcitons in quantum wells 

The bound states between two charged fermions, an electron from the conductive band and a 
hole from the valence band, in the presence of a magnetic field are called magnetoexcitons. 
In what follows we consider a single quantum well (SQW) and coupled quantum wells 
(CQW's) made with direct-gap semiconductor that has nondegenerate and isotropic bands: 
Ec{k,kz) = Eg + h\^ /2mc + h^k'^/2mcand Ey(k,k;,) = /2my + h^k'^/2my, where 
k is a two-dimensional (2D) wave vector. Eg is the semiconductor band gap, and rric (m^) 
is the electron (hole) effective mass. The z-axis is chosen to be the axis of growth of the 
quantum-well structure, and the constant magnetic fields is B = (0, 0, B). The x-y plane 
has been taken to be the plane of confinement. In what follows we neglect any electron-hole 
correlations along the z-axis. This approximation takes place when the effective mass of the 
hole considerably exceeds that of the electron and the slow motion of the hole is separated 
from the fast motion of the electron. The assumption is applicable for many crystals of 
A^^^B^ type. In the presence of confinement potentials Uc,v{z), the corresponding electron 
ip and hole (j) wave functions are defined by the solutions of the one-particle Schrodinger 
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equations: 

Here, Exc{E^y) is the electron (or hole) confinement energy, A and ^ denote the quantum 
numbers of the states in the confinement potential. For simplicity, we shall take into ac- 
count only the first electron Eqc and hole Eq^ confinement levels. In the above equation Zc 
and are the electron and hole z-coordinates, respectively. 

The exciton motion in (x, y)-plane changes its spectrum, i.e. the magnetoexciton energy 
E{Q) depends on the in-plane exciton pseudomomentum ?iQ = h{Qx, Qy,0)- The influ- 
ence of the exciton motion on its spectrum (the dispersion relation), in the case of a SQW 
or CQW's have been studied extensively over the past decades [1-5], assuming that the 
magnetoexciton energy -E'(Q) and the corresponding wave functions of the relative internal 
motion can be obtained from the solutions of the Schrodinger equation with the following 
Hamiltonian: 



H = Eg + Eo, + Eo, - —Vl + B X r .Vr + - Vc{r + Rq • (1) 
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Here, r = Tc — r^,, and r^j-^,) is the two-dimensional (2D) electron (hole) position vector, 
/i = mcniy/M is the exciton reduced mass, where M = nic + is the exciton in-plane 
mass. 7 = {m^ — mc)/M, Rq = ^^Qo> where Qg = (— Qy,Qa;,0), and I = (JicjeBYl'^ 
is the magnetic length. Vq represents the electron-hole Coulomb attraction screened by the 
high-frequency dielectric constant Coo : 

where the structure factor / is defined by: 

/ + 00 /■ + 00 

dzcl dz^Qi^^{-q{zc- z^y\\(f\j,z^<\)l^{z^). (3) 
-OO J —CO 

Since the last term in ([T|) is the only term which depends on the exciton momentum, the 
following statements take place: 

(i) The magnetoexciton dispersion does not depend on the electron and hole masses. 

(ii) The magnetoexciton dispersion is determined only by Coulomb interaction. 



2 Dimensional reduction in the dynamics of bulk magnetoexci- 
tons 

Strictly speaking, the excitons are bound states between two charged fermions, and there- 
fore, the appropriate framework for the description of the bound states is the Bethe-Salpeter 
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(B-S) formalism [6-9]. In the absence of a magnetic field, by using a series of approxi- 
mations (such as the introduction of the equal-time wave function, the assumption that the 
B-S kernel depends only on the difference of the relative momenta) the B-S equation for 
electron-hole bound states can be simplified to the well-known Schrodinger equation for 
the relative internal motion [10]. The existence of a magnetic field induces a coupling be- 
tween the center-of-mass and the relative internal motions, because even a small transverse 
exciton velocity (or small transverse wave vector Q) will induce an electric field in the rest 
frame of the exciton which will push the electron and the hole apart, so the binding energy 
must decrease as the transverse velocity increases. Thus, one can expect that in the presence 
of a magnetic field the simplification of the B-S equation to the Schrodingier equation is 
not trivial. 

Several non-trivial effects produced by magnetic fields have been recently predicted in 
quantum field theories. For example, in the massless QED, by means of the lowest Landau 
level (LLL) approximation, the B-S equation has been reduced to the Schrodingier equa- 
tion, and as a result, it was predicted that the external constant magnetic field generates an 
energy gap (dynamical mass) in the spectrum of massless fermions for any arbitrary weak 
attractive interaction between fermions [11-14]. This effect is model independent (univer- 
sal), because the physical reason of this effect lies in the fact that dynamics of the LLL is 
essentially D — 2-dimensional. In other words, the essence of this effect is the dimensional 
reduction (from 3 + 1 to 1 + 1, or 2 + 1 — > + 1) in the dynamics of fermion pairing in the 
presence of a the constant magnetic field. Later, it was suggested that a similar effect could 
explain some experimental findings in the physics of high-temperature layered supercon- 
ductors [15]. In what follows, we will see that the dimensional reduction in the dynamics 
of magnetoexcitons manifests itself in the fact that the magnetoexciton dispersion does not 
depend on the electron and hole masses. 

We first use the B-S formalism to describe excitons in a bulk material in the presence of a 
strong constant magnetic field B along the z-axis. After that, we apply the bulk B-S for- 
malism to a SQW or CQW's. The process of generalizing the bulk equations to the case 
of quantum-well structures is a straightforward procedure because of the assumption that 
there are no electron-hole correlations along the z-axis. 

The basic assumption in the B-S formalism is that the electron-hole bound states are de- 
scribed by the B-S wave function (B-S amplitude) 2) = *(rc, r^; Zc, Zv;ti,t2), where 
the variables 1 and 2 represent the corresponding coordinates and the time variables. This 
function determines the probability amplitude to find the electron at the point (rc, Zc) at the 
moment ti and the hole at the point (r^, z^) at the moment t2. The B-S amplitude satisfies 
the following equation: 



Here / is the irreducible B-S kernel, and Gc,v are the electron and the hole Green's func- 
tions. If the screening effects are taken into account by means of the high-frequency dielec- 



*(1;2) 
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(4) 
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trie constant 



£00, then the irreducible kernel is given by 



T// ^ ^^^^ / 

V{r;z) = -—- / 

too J 



(27r)2 2tt |q|2 + g2 



d^q dqz 1 



exp (q.r + q^z)] . 



(5) 



In what follows, we use the center-of-mass (R, Z) = (acfc + ctyVy, UcZc + a^Zy) and 
the relative (r, z) = (rc — rv,Zc — Zy) coordinates. The coefficients etc = (1 — 7)/2, 
at, = (1 + 7)/2 are expressed in terms of the parameter 7 = (m-u — rric) / {nic + mt,,) which 
accounts for the difference between the electron and the hole masses. The B-S equation for 
the equal-time B-S amphtude in the center-of-mass and reduced coordinates assumes the 
form: 



The B-S amplitude depends on the relative internal time t — t' and on the "center-of-mass" 
time: 



^q q^ (r, R; z, Z- 1, t') = exp ( ^|^(act + a^t') J ^-q^q^ (r, R; z, Z; t-t'), (7) 



where E{Q, Qz) is the exciton dispersion. Introducing the time Fourier-transforms accord- 
ing to the rule f{t) = f{uj) exp {icjt) we transform the B-S equation into the 
following form: 



ten as a product of phase factors and translation invariant parts. The phase factor depends 
on the gauge. In the symmetric gauge the vector potential of the magnetic field A is defined 
by A(r) = (1/2)B x r, and the Green's functions are [16]: 



^QQ^{r,R;z,Z;t,t) = f dz' dZ' cfr' cfR' dtidt2 
Gc{R + a„r, R' + Q,r'; Z + ^z, Z' + ^z'; t - h) 
G„(R' - ay, R - aer; Z' - Z - ^z; h - t) 



V{r'-z')^QQ^ir',R';z',Z';h,h). 



(6) 




^Q^Q^{r,R- z,Z;lo) = J dz' dZ' dFr' d^R! ^ 
Gc (R + ayV, R' + a^; Z + a-^z, Z' + a^z'; Uu + acE{Q, Q^)) 
Gy {R! — act' , R — act; Z' — acz' , Z — acz; hu — avE{Q, Q^)) 




Gc,„(r,r';2;, 2';a;) = exp r.A(r') Gc,v{r -r';z - z';u). 



(9) 



The broken translation symmetry requires a phase factor for the B-S amplitude: 



f^)=exp z— r.A(R) XQ,Q,(r,R; 2, ^; ^^)- 



(10) 
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The B-S equation ([8]) admits translation invariant solution of the form: 
The function xq q (r; ^; satisfies the following B-S equation: 



(11) 



XQ^Q^ (r; z;u;)=J dz'dZ'd^r'd^R'^ exp [f ((r + r').A(R' - R) + 7r.A(r'))] 

Gc(R - R' + a^,(r - r'); Z - Z' + ay{z - z'); huj + OcE) 
Gy(R' - R + adr - r'); Z' - Z + a^z - z'); huo - a^E) 

y(r';z')XQ,Qjr';z';17). (12) 

The substitution R' — > R' + R + 7r provides the following equation 
for the Fourier transform of the exciton wave function XQ g (^]kz]uj) = 
J dzd'^rex.p —i (k.r + kzz) xq q (r; z] ^) of the exciton wave function: 

m,QS^ - ^) = / ^ ^ exp [-^{q + Q).R] x 

Gc (iq + k - ^A(R); + a^Qz] Tiuj + Oc^) x 
Gy (-iq + k - ^A(R); kz - UcQz', fioJ - a^E^ x 
f(p- [k-lfA(R)] ■,Pz-kz)xQ,Q^{v-lQ:Pz;^), (13) 

where y(k; /c^) = — (47re^/eoo) (k^ + k1) ^ and Gc,„ (k; A;^; Titj) ai-e the Fourier trans- 
forms of Gc,v (r; z; fiu). 

In the effective-mass approximation the exact fermion Green's functions Gc^v are replaced 
by the corresponding propagator of the free fermions G'f}, . The translation invariant parts 
G^}, can be decomposed over the Landau level poles: 

Gi^ (r; z- huj) = ^^d^Gf), (k; fc,; hoj) exp i (k.r + kzz) , 



GS(k; kz-, huj) = 2 E^=o(-l)" exp (-/^k^) L„ (2/% 



n^kl/2mc + Eg + nVLc{n + 1/2) 



-1 



Gl°)(k;A;,;nL^) = 2E^=o(-irexp(-/2k2)L„ (2/2k2) x 



(?la;+ 7i^fc2/2m^, + + 1/2) 



20+ 



-1 



(14) 



Here Ln{x) are the Laguerre polynomials, and fi^c,v = heB/cmc,v are the electron and 
hole cyclotron energies. In strong magnetic fields the probability for transitions to the 
excited Landau levels due to the Coulomb interaction is small. Thus, the contributions to 
the Green's functions from the excited Landau levels is negligible, and therefore, one can 
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apply the lowest Landau level (LLL) approximation, where we keep only n = term in 



G^{k;k^;huj) 2 exp (-/V) (hu 



Eg + n^ki/2mc + hQ.c/2 



+ 



-1 



(15) 



The solution of the B-S equation in the LLL approximation can be written in the following 
form: 



iRo.k 



(16) 



Thus, the LLL approximation reduces the problem from 3 + 1 dimensions to 1 + 1 di- 
mensions problem for obtaining functions <^Q^{kz;u}) and the energy E{Q, Qz) from the 
following equation: 



nuj+acE-{Eg+^{k:,+acQz)^ + 



+ 



.(17) 



In the LLL approximation, the in-plane exciton dispersion is determined by the Coulomb 
interaction: 



d^q 2 / .exp[zq.(r + Ro)] 



rV'oo(r) 



(18) 



Soo J (2^)^'^'"' ' {q^ + kl) ■ 
Here, V'oo(r) = exp (— is the ground-state wave function of a hydrogen atom 
in a magnetic field. The solution of ([TT] ) can be chosen in the following form: 



huj + OcE 



-1 



hto — a^E + 



^{kz-a.Qz? + ^ 



^0^ 



(19) 



where cpg^ (kz) is a function to be determined. By integrating both sides of ([TSl l over u, we 
find the following equation for the exciton wave function 

^Q.ikz) = I ^^Q^{kz,u;) = cPqM/ {E-Eg- h'kl/2f, - T^QI/2M) 

and exciton energy Eh{Q,Qz) = Eg + }-hQ. — E{Q,Qz) = heB / fi is the exciton 
cyclotron energy): 

+ E,{Q,Qz)UQAkz)- J -^lQ{kz-Pz)^QApz). (20) 







2/x 



2M 



The exciton binding energy Ef, > could be obtained from the solutions of (120]) by means 
ofEh = Eb{Q = 0,Qz = 0). 

In the case when Q = and Qz = 0, eq. (l20l ) is similar to the well-known one-dimensional 
Schrodinger equation for a hydrogen atom in the adiabatic approximation [17-20]. 
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3 Magnetoexciton dispersion in quantum wells in the lowest 
Landau level approximation 

The assumptions that: (i) we neglect any electron-hole correlations along the z-axis, and (ii) 
we take into account only the first electron Eq^ and hole Eq^ confinement levels with wave 
functions (^oc(-Zc) and (j)Qy{zv), respectively, greatly simplify the description of the motion 
along the z-axis. In the cases of a SQW and CQW's, the Fourier transform of the exciton 
wave function satisfies the following B-S equation: 

XQ(k - iQ;u;) = J ^^d^Rj^^ f exp [-^{q + Q).R] x 



Gc (h + - ^A(R); no; + o^e) (-iq + k - ^A(R); huj - a,E) 



2e 

he 



A(R) 



XQ(p-iQ;J7), 



(21) 



where the potential l^(k) = — (27re^/(|k|)/eoo) |k|^^ depends on the quantum-well ge- 
ometry through the structure factor/ (k). 

In the LLL approximation the exact fermion Green's functions Gc,i> are replaced by the 
corresponding propagator of the free fermions G^c)j '■ 



Gc(k; huj) ^ 2 exp {-Pk^) {Tito - [Eg + Eq^ + n9.c/2] + iO+) ^ , 
G„(k; huj) « 2 exp (-/^k^) {huj + Eq^ + - «0+)"^ . 



(22) 



The solution of the B-S equation in the LLL approximation can be written in the following 
form: 



'k:uj) 



exp 



zRo.k 



(23) 



Thus, the LLL approximation reduces the problem from 2 + 1 -dimensions to + 1- 
dimension problem. The function energy E{Q) can be obtained from the following 

B-S equation: 

{huj + ttcE -Eg- Eoc - nnc/2 + lO+y^ {hoj - a^E + Eov + - t0+y\24) 

In the LLL approximation, the exciton dispersion is determined by the term: 



m) 



27re2 



d\ 



(2vr)' 



/(|q|)exp bq.(r + Ro)] 



(25) 



The solution ^e{^) of (l24l) can be chosen in the following form: 



huj + acE - Eg - Eqc 



2 



+ 20+) (hio - a^E + Eoy + 



-1 



• (26) 
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Integrating both sides of B-S equation (l24l) over lo, we find that the exciton dispersion is 
determined only by the Coulomb interaction (1251 ): 

^(IQI) =Eg + Eo, + Eo, + nn/2 - /(IQI). (27) 

It turns out that in the LLL approximation the magnetoexciton dispersion does not depend 
on the electron and hole masses and is determined only by Coulomb interaction. 
The LLL approximation greatly simplifies the equations, but we may ask whether the mag- 
netoexciton dispersion will be significantly affected by the contributions from the infinity 
number of Landau levels with indexes n > 1 neglected in the LLL approximation. In the 
next Section we address this question. 



4 Magnetoexciton dispersion in GaAs — {Ga, Al)As quantum 
wells 

In the previous two Sections, we decomposed the single-particle electron (hole) Green's 
function over the Landau poles and we kept only the term with index n = 0. This term is 
relatively simple, and allows us to perform all integrations in the B-S equation ([T3] ). Unfor- 
tunately, the terms with n > 1 are more complicated, and it is impossible to perform the 
integrations over the corresponding variables. 

There exists another approach which allows us to figure out the contributions to magne- 
toexciton dispersion due to the Landau levels with indexes n > 1. It starts from the B-S 
equation but rewritten in the following form [21, 22]: 

{^^mi-E, - 2k [-^^^r. + fA(x„ye,.e)]' - - U,{zc)) x 

= iV{rc-r^;Zc - Zt,)^'(rc, r^,; z^, z^; ti, ti), 

where V{r, z) is defined by Since there are no electron-hole correlations along the 
z-axis, we separate the variables and write the B-S amplitude in the following form: 

*(rc,2:c,ii;r^,,2:„,t2) = exp|z Q.R - ^r.A(R) - f (ac^i + at.i2) }x 

XQ{T^;ti-'t2)^Q{zc)(t}o{zv), (28) 

where E = E{Q) is the magnetoexciton dispersion. After some tedious, but straightfor- 
ward calculations, we arrive at the conclusion that the Fourier transform of the B-S ampli- 
tude 



/(pn r+00 ^ 



(29) 
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satisfies tlie following equation [21,22]: 



J§$Jdh exp (z(q' - q).r) \hn - n,{q', Q) - nf^ (Q, q'; r) 



nn - n,{q', Q) - (Q, q'; r)J xqCq'; 
. . d2q' 2^eV(|q-q'|) ^f^' v^^n'- O'^ 



(30) 



Here, we use the following notations: 
J7e(q,Q) = Ec{q+acQ) + Eoc-acE, n^{q,Q) = -E^{q-a^Q)-Eo^+ayE, (31) 

eh , „ eh ,„ , e^i?^ 



J^?MQ,q;r) 



2Mc 



B± X r).Q + (Bx X r).q + ^H, 



r!^(Q,q;r) = ;^(Bxxr).Q 



2mcC 
eh 



;B± X r).q + 



(32) 



(33) 



where Ec^y{q) = Ec^^{q, = 0). We are looking for the solution of Eq. ( [30l ) of the form: 



XQ(q;f^) 



5Q(q) 



[/i^? - f]c(q, Q) + iO+] [nf] - r?„(q, Q) - iO+] ' 



(34) 



where 5'Q(q) is a function to be determined. 

We introduce the function XQ(q)» which is the Fourier transform of the equal-time B-S 
amplitude (or exciton wave function) XQ(r) = XQ(r; ti — ^2 = 0): 



XQ(q) 



-oo 2tt 



XQ(q;f^) 



(35) 



By taking into account the analytic properties of xq{(Ii^)^ we obtain the following B-S 
equation for determining the exciton energy E' = E{Q) — Eg — Eqc — Eq^ and the Fourier 
transform of the exciton wave function XQ(q)- 



/ iSji [(# + ^3!^) - 1') + "?(Q. 1. 1') + n?(Q.<i. 1') - '-^^ 

In what follows, the last term in 061 ) will be referred as the B-S term: 



(36) 



Vr. c(n n'-n r'\ - [EvW-cvQ)-E.{q-avQ)]n^{Q,q,q') 
VB^siq, q , ^ ) - E'-E.iq'+a^QyE^q^a.Q) 



+ 



[E,{q'+a,Q)-E,{q+a,Q)]n^{Q,q,q') 

E'~Ec{q+arQ)~E„{q'~a^Q) 



+0B(Q,q,q') 



E'-Ec{q'+acQ)-E^{q-a^Q) E'-Ec{q+acQ)-E^{q'-ayQ) 



■ (37) 



10 



Here, the following notations have been used: 

f^?.(Q,q,q') = / ^^'rexp [z(q' - q).r] (Q^q/.^), (38) 

^^?.(Q,q,q') = / d'rexp[<q'-q).r] f]?(Q,q';r)l}?(Q,q';r). (39) 

In position representation, the B-S term generates a non-local potential which depends on 
the energy E': 

VB-s{ry;Q,E') = J ^ J ^yB_s(q,q';Q,i?')exp[z(q.r-q'.r')]. (40) 

The solution of Eq. 061 ) can be written as 

Xqiq) = exp (-iq.Ro) (q - Qo) , 
where the function ^ (q) satisfies the following equation: 

E'^{q) = (q) - ^g (B^ x q) .Vq* (q) - ^Vl^ (q) 

~~ J WF P inq - q J -koj |q_q'| ^iq J 

- / exp [t (q - q') .Rq] Ve^siq + iQo, q' + iQo; Q, i^')^(q')- (41) 

The B-S equation (|4T]) differs from the Schrodinger equation. If we neglect the B-S term 
in the right-hand side of (|4T]) . we obtain the Schrodinger equation for magnetoexcitons 
with the Hamiltonian ([T]). It can be seen that according to the Schrodinger equation, the 
magnetoexciton dispersion is totally determined by the Coulomb term, while according 
to the B-S equation, the effective potential (l40l ) also contributes to the magnetoexciton 
dispersion. 

Since the Bethe-Salpeter term plays an important role in determining the magnetoexciton 
dispersion (see the next two Sections), one may well ask a question about the physical 
meaning of this term. The answer is that the B-S term takes into account the contributions 
to the single-particle Green's functions (fT4l) from the Landau levels with n > 1. 

5 Magnetoexciton dispersion in single GaAs/AlxGai-xAs 
quantum well 

In this Section, we first calculate the ground-state energy of a heavy-hole magnetoexciton 
with a zero wave vector (Q = 0), assuming a single GaAs quantum well with a thickness 
L sandwiched between two Al^Gai-xAs layers. The electron in-plane mass rric and the 
electron z-mass rricz are chosen to be rric = rricz = 0.067mo, where niQ is the bare electron 
mass. The in-plane heavy-hole mass m^, and the hole z-mass rrivz are expressed in terms 
of the Luttinger parameters 71 and 72: m„ = mo/ (71 + 72) and m^^ = mo/ (71 — 272). 



11 



L{nm) 


B{T) 


/5 


Eco{meV) 


E^o{meV) 




EexpieV) 


EsieV) 


4.03 





0.786 


100 


26.9 


1.6355 


1.638 


1.6355 


4.03 


2 


0.810 


100 


26.9 


1.6356 


1.639 


1.6357 


4.03 


4 


0.869 


100 


26.9 


1.6365 


1.640 


1.6367 


4.32 





0.776 


93.5 


24.3 


1.6262 


1.630 


1.6262 


4.32 


2 


0.802 


93.5 


24.3 


1.6265 


1.631 


1.6266 


4.32 


4 


0.861 


93.5 


24.3 


1.6274 


1.632 


1.6275 


7.2 





0.702 


51.0 


11.0 


1.5716 


1.571 


1.5716 


7.2 


2 


0.734 


51.0 


11.0 


1.5719 


1.572 


1.5720 


7.2 


4 


0.803 


51.0 


11.0 


1.5730 


1.573 


1.5731 



Table 1 : Variational calculations of the heavy-hole exciton ground-state energies with Q = 
for various well widths L and weak magnetic fields B. The trial function (l42l) depends 
on the variational parameter /?. The energy gap is Eg = 1.519 eV. The electron and hole 
confinement energy levels Eco and E,uq are calculated assuming squared-well potentials of 
finite depths. The E'^ja^-column represents the results from the variational calculations with 
the following Luttinger parameters: 71 = 7.36 and 72 = 2.57 [23]. The measured ground 
state energies E^xp are reproduced from [24]. The £^s-column represents the ground-state 
energies calculated according to the Schrodinger equation with the Hamiltonian ([D 



It is known that the difference between the bandgap energies of GaAs and AlxGai-xAs 
provides a finite potential well, confining the electron-hole pairs in the Galas quantum well. 
We assume that the potentials are square-well potentials of finite depths Vc = 0.6AEg{x) 
and Vy = OAAEg{x), respectively. The energy-band-gap discontinuity [23] is assumed to 
be AEg{x) = (1.555x + 0.37x'^)meV. The confinement energy levels Eco and E^o are 
obtained by solving the following transcendental equations: 



tan 



tan , 

Here, Eb = Tl^ /2^a\ is the exciton Bohr energy. The structure factor f{k) is calculated 
by means of the following wave functions: 




as 



^cA^) = B,,, COS ( • / ^-■-^^o,.o} 



Ij,Eb 



-00 < Z < -1/2, 
-1/2 <z< 1/2, 
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c,v 



Brv exp 



L / 'm-cz,vz{Vc,v — EcO,vo) 

2aB V IJ-^B 



COS 



L nT-cz,vzEcO,vQ) 

2aB V M^^s 



(42) 



Since the B-S equation (|4TI) is rather complicated we shall obtain numerical results for the 
ground-state energy within the framework of the variational approach. In the case of weak 
magnetic fields, i.e. << Eb, we use a hydrogen-like trial function with a variational 
parameter (3: 

2^2(3 ( 2r(5\ 

With this trial function we calculate the following magnetoexciton energy: 

E = Eg + Eco + EyQ — E{P)Eb, 
where E{P) is defined by the solution of the following equation: 

3 / ^ 2 



E{P) = -4/?2 + l28/?3 
With the trial function 



dx- 



(16/32 + 2;2)3/2 128/32 \Eb 
, the B-S contribution to the ground state is 



+Vb-s{P,E,B), (43) 



+8a|E4/36(-195 + 896-^2 - 3672) - Aa^^E^ (3^ {-17?, + 128^2 ^ 4^2) 
+64a|5£;2^i0(4i _^ 1408£;2 _ 322^2 _ 492^4^ 

-512a|S/3i2(3 + 208^2 _ ig^2 ^ ;l5^4) _ 32a|E3^8(79 + 1152^2 ^ §0272 
+1727^) + 1024/3i4[48^2 ^ ^ ^2)2]] 
-64^2^8 j_2048a|£;/36 + 1024/38 + 4&a%l3^{l + 32^2 + ^2 _ ^2^4) 
-lQa%Ef3'^{'i + 32£;2 _^ 2472(2 + 72) 



-Fa|^2(g4£;2 _ 3[;li ^ 8^2(7 ^ ^2)]^] 



alE 



}• 



(44) 



The dimensionless variables E and in the right-hand side of Eq.(l44l) must be replaced 
by E{f3)hQ./ E1^ and ub/I, respectively. The results obtained by using the hydrogen-like 
trial function are presented in Table 1 . We used more significant figures to stress on the fact 
that the magnetoexciton energies, calculated by applying the B-S formalism are extremely 
closed to those, provided by the Schrodinger equation. 

The magnetoexciton dispersion are determined by the Coulomb interaction and the B-S 
term in Eq. (|4TI ). The contribution from the Coulomb interaction to the energy of the mag- 
netoexciton (in Eb units) increases quadratically for small wave vectors QaB << 1, and 
can be written as {QaB)'^fJ-/Mc- The hydrogen-like trial function provides the following 
expression for the in-plane exciton mass Mc- 



Mc 



32/3^ 



RV 



dx 



x2/(x-^) 



ubJ Jo (I6/32 +x2)3/2' 
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L{nm) 


B{T) 


/3 




EexpieV) 


EsieV) 


Mc/mo 


MB-s/rr. 


4.03 


20 


0.85 


1.650 


1.644 


1.651 


0.145 


0.0025 


4.03 


18 


0.84 


1.648 


1.643 


1.649 


0.127 


0.0010 


4.03 


16 


0.84 


1.647 


1.642 


1.647 


0.114 


0.0002 


4.32 


20 


0.84 


1.641 


1.636 


1.642 


0.147 


0.0026 


4.32 


18 


0.83 


1.639 


1.635 


1.640 


0.129 


0.0011 


4.32 


16 


0.83 


1.638 


1.634 


1.638 


0.116 


0.0002 


7.2 


20 


0.86 


1.587 


1.583 


1.588 


0.176 


0.0044 


7.2 


18 


0.84 


1.585 


1.582 


1.586 


0.159 


0.0022 


7.2 


16 


0.84 


1.583 


1.581 


1.584 


0.142 


0.0007 


7.49 


20 


0.86 


1.584 


1.580 


1.584 


0.178 


0.0046 


7.49 


18 


0.84 


1.582 


1.579 


1.582 


0.161 


0.0024 


7.49 


16 


0.84 


1.580 


1.578 


1.580 


0.144 


0.0008 


7.5 


14.5 


0.67 


1.577 


1.577 


1.572 


0.131 


0.0302 


7.5 


12 


0.64 


1.575 


1.573 


1.570 


0.049 


0.0160 


7.5 


8.5 


0.60 


1.572 


1.570 


1.569 


0.026 


0.0071 



Table 2: Variational calculations of the heavy-hole exciton ground-state energies for var- 
ious well widths L and strong magnetic fields B. The trial function (1451 ) depends on the 
variational parameter (3. The energy gap is Eg = 1.519 eV for the L = 4.03,4.32,7.2, 
and 7.49-nm wells, and Eg = 1.512 eV for the L = 7.5-nm. The Eyar-colurrm repre- 
sents the energies obtained by the variational method using the following Luttinger param- 
eters: 7i = 6.9 and 72 = 2.4 [25]. The measured ground state energies E^xp for the 
L = 4.03, 4.32, 7.2, and 7.49-nm wells are reproduced from [24], and for the L = 7.5-nm 
well from [26]. The £'5-column represents the ground-state energies calculated according 
to the Schrodinger equation. The Mc and Mb^s are the masses calculated according to 
Eqs. ^ and 
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The contribution to the exciton dispersion due to the B-S term can be evaluated analytically. 
We found that it also increases quadratically for small wave vectors, but for B < AT, this 
contribution is about one tenth of {QaB)'^fJ'/Mc- Thus, in a weak magnetic field, there is no 
measurable difference between the results calculated by the Schrodinger equation, and these 
obtained by the more complicated B-S formalism. For a weak perpendicular magnetic field 
and small wave vectors, the Coulomb interaction dominates, which means that a hydrogen 
type of ground state slightly modified by the magnetic field exists. 

Next, we consider the case of a strong magnetic field. In this regime we choose the trial 
wave function ip/sir) to be similar to the corresponding ground-state wave function of a 
charge particle in a magnetic field, but depending on a variational parameter (3: 



Here, and in what follows, we use the exciton cyclotron energy hO, for energy unit and mag- 
netic length R for unit length. The ground state magnetoexciton energy will be calculated 
by minimizing the energy functional E'{P) = {E — Eg — Eq^ — EQ^)/h^ with respect to 
the variational parameter (3: 



Note, that (i) all terms in the last equation are dimensionless (in a cyclotron energy HQ 
unit), and (ii) we have written the contributions from the Coulomb interaction and from the 
B-S term ( [37] ) as a sum of Q-independent terms, Vc{P) and Vb-s{P, E'), and Q-dependent 
terms, Vc{P,Q) and Vb-s{/3, E' ,Q). The Q-dependent terms will be used to obtain the 
magnetoexciton dispersion. The second and the third term in ( |46l) are given by: 



VB-siP, E') = SZ^^^^-l±2^{e^E'P' [-56E'^P'j^ + 32E'^P<^j^ + (-1 + 7^)^ 

-32£;'2 [-1 + /34^2 ^ (3 ^ AE'P^{-2 + E'P^))-/'^]] Ei(4^'/32)}. (48) 

Here, E^ = ^Ji^^e^ j (eoo-R) is the binding energy of the two-dimensional (L = 0, /3 = 1) 
magnetoexciton, calculated according to the Schrodinger equation. 

The energy of the magnetoexciton increases quadratically for small wave vectors 
{QR « 1): Vc{P,Q) = [fi/2Mc{L,B,P)]{QR)^ and Vb-s{P,E',Q) = 
[ii/2Mb-s{L, B, P)]{QRf. The in-plane mass Mc{L, B, P) is due to the Coulomb inter- 
action and does not depend on the electron or the hole mass: 




(45) 





(47) 




(49) 



15 



where M2D = 2^/^eoo?i^/ {y/'Ke^R). The second in-plane mass, Mb-s, has its origin in the 
fact that the B-S term depends on Q, and for QR << 1, M b-s is defined by the following 
equation: 

Wifim) = ^^^^Sfe^{e^^'^^[256i?'5/3i2/ + 64i?'V^°7^(5 - 177^) 

-3/?2(-l + 72)3 _ 2E'P^{-1 + 72)2(1 + 1272) - 48S'2/?6^2(2 _ 7^2 ^ 5^4) 

+lQE'^[-2 + 2/34 + /?S^(4 - 5372 + 747^)]] 

-64£;'3/32[_/j2 ^ ^ 4^/2^10^4(5 _ 18^2) ^ ^6^2(_7 ^ 33^2 _ 30^4) 

+2E' [-1 + + ^S^(2 - 2972 + 457^)]]Ei(4^'^2)|^ (50) 

where Ei(2;) = — dt exp(— is the exponential integral function (the principle value 
of the integral is taken). 

Table 2 gives the results of our variational calculations. It can be seen that the B-S equation 
provides similar results for the ground-state energies as the Schrodinger equation does. 
Since the B-S mass is much smaller than the Coulomb mass, one can say that in strong 
magnetic fields the exciton dispersion for small wave vectors {QR « 1) is determined by 
the B-S term rather than the Coulomb interaction. 



6 Coupled quantum wells in strong magnetic fields 

In this Section, we consider exactly the same double well electron-hole system as in Refs. 
[27,28]. The electron layer and hole layer have finite widths, denoted below by Lc and L^,, 
and they are separated by a distance D. We assume that the electrons and holes are confined 
between two parallel, infinitely high potential barriers. This assumption greatly simpli- 
fies our numerical calculations of the magnetoexciton energy and the Coulomb mass, but 
by neglecting the existence of the finite confinement potentials, we cannot provide a more 
realistic value for this part of the exciton energy related to the exciton confinement along z- 
direction, than the sum of the well-known terms fi^i:'^ /2mc^vLl ^. Obviously, the more real- 
istic model of a symmetric (or asymmetric) DQW with finite quantum-well widths [29, 30] 
will cause minor corrections to our main conclusions, which are: (1) the B-S formalism 
provides a term, which does not exists in the Schrodinger equation, and (2) the term plays 
an important role in determining the magnetoexciton dispersion. 

The basic features of the CQW's magnetoexcitons are the same as that of the SQW mag- 
netoexcitons. However, because of the separation between the electron and hole layers, 
the Coulomb energy and the Coulomb in-plane mass differ quantitatively from those of the 
SQW magnetoexciton. In other words, in strong magnetic fields, Eq. (l46l ) holds, but the 
Coulomb interaction and the corresponding in-plane mass are defined as follows: 



= V?/o°° dxe-^-^F (x, ^, ^) , (51) 
M2D _ ^/I foo , 2 p(x ^] (52) 



16 



B{T) (3 E,ar{meV) Es{meV) ^ 



Mb-S 



10 0.96 6.36 6.56 2.06 

9 0.96 5.17 5.43 1.75 

8 0.96 4.03 4.31 1.46 

7 0.96 2.94 3.20 1.19 

6 0.95 1.91 2.11 0.95 

5 0.94 0.95 1.04 0.72 

4 0.92 0.01 0.01 0.52 



0.228 
0.221 
0.216 
0.215 
0.218 
0.230 
0.247 



Table 3: Variational calculations of the magnetoexciton energies for various strong mag- 
netic fields B, measured relatively to the Eg + Eqc + Eq^ level. The trial function (1451 ) 
depends on the variational parameter /3. The E^^q,. -column contains the energies calculated 
by the variational method with the following parameters: rric = 0.067mo, m^, = O.lSmo, 
eoo = 12.35, Lc = Ly = 8nm, D = 11.5nm. The E^g-column represents the magnetoex- 
citon energies calculated according to the Schrodinger equation. Mc is the in-plane mass 
defined by Eq. (l52l) . The Mb-s is the mass calculated according to Eq. (l50l) . 

In CQW's, the structure factor is: 



Table 3 gives the result of our numerical calculation of the magnetoexciton energy, but rel- 
atively to the Eg + Eqc + Eqv level. We used the same parameters as in Refs. [27] and [31]. 
It can be seen that the B-S equation provides slightly different results for the binding energy 
than the Schrodingier equation. 

The main difference between the B-S and the Schrodinger equation is in their predic- 
tions about the in-plane magnetoexciton mass in a strong magnetic field. Unfortunately, 
optical experimental studies can provide information about the exciton dispersion only for 
Q < Qph, where hQph is the photon momentum. Other studies, such as the photolumines- 
cence measurement experiments which can measure the exciton-mass dependence of the re- 
combination time, or experimental data related to the polariton effects, can provide informa- 
tion about the magnetoexciton dispersion. Many of these experimental techniques [32-36] 
are used to measure the magnetoexciton dispersion in the presence of an in-plane magnetic 
field. As we mentioned above, the measurable differences between the magnetoexciton 
dispersions, as predicted by the B-S formalism and by the Schrodinger equation, are to be 
expected in strong perpendicular magnetic fields. To the best of our knowledge, there is 
only one paper [27] where the exciton dispersion in Ga^s/Gao. 67^^0.33^-5 CQW's in a 
weak perpendicular magnetic field has been measured. There is a good agreement between 
the mass Mc and the measured mass in a weak magnetic field. Referring to the conclusion 
that the B-S term in a weak magnetic field has a very small contribution to the dispersion 
compare to the contribution due to the Coulomb interaction, one can say that there exists a 
good agreement between the B-S formalism and the measurements. 



1671^(1 - e-«=^)(l - e-?''^)e-'^^ 
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Next, we discuss the fact that Mc increases by about 4 times if we increase the magnetic 
field from 4T to lOT. If the magnetoexciton dispersion in strong magnetic fields (B > 5T) 
is determined mainly by the B-S term, then the magnetoexciton mass should not increase 
so dramatically, and therefore, new experimental points are needed to prove or disprove the 
conclusions drawn by applying the B-S formalism. 

7 Conclusion 

We have applied the B-S formalism to the quantum-well excitons in a constant magnetic 
field applied along the axis of growth of the quantum-well structure. We found that (1) 
in the LLL approximation the B-S equations provides the same results as the Schrodinger 
equation; (2) beyond the LLL approximation, the B-S equation contains an extra term (B-S 
term). This term takes into account the transitions to the Landau levels with indexes n > 1. 
We applied a variational procedure to obtain the effect of the B-S term on the magnetoex- 
citon ground-state energy and magnetoexciton mass. We used a simple hydrogen-Uke trial 
wave function in a weak magnetic field, and figured out that in a weak perpendicular mag- 
netic field the results obtained by the B-S formalism are very close to the results calculated 
by means of the Schrodinger equation. In a strong magnetic field, we used a trial function 
similar to the wave function of a charged particle in a magnetic field. We calculated that in 
a strong magnetic field, the ground-state energy is very close to that obtained by means of 
the Schrodinger equation, but the magnetoexciton dispersion is determined by the B-S term 
rather than the electron-hole Coulomb term in the Schrodinger equation. 
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